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Abstract 

Light clusters (mass number A < 4) in nuclear matter at subsaturation 
densities are described using a quantum statistical approach. In addition 
to self-energy and Pauli-blocking, effects of continuum correlations are taken 
into account to calculate the quasiparticle properties and abundances of light 
elements. Medium-modihed quasiparticle properties are important ingredi¬ 
ents to derive a nuclear matter equation of state applicable in the entire 
region of warm dense matter below saturation density. The influence of the 
nucleon-nucleon interaction on the quasiparticle shift is discussed. 


1 Introduction 

The nuclear matter equation of state (EOS) is an important ingredient to 
understand the properties of nuclear systems, for instance for heavy ion colli¬ 
sions (HIC), see [T], and for astrophysical applications such as the formation 
of neutron stars in core-collapse supernovae, see [2]. In both cases, matter is 
considered at subsaturation density (baryon density ns < n-sat ~ 0.16 fm“^), 
moderate temperatures T < 20 MeV, and proton fraction Yp = /ub be¬ 
tween 0 and 1 where n*”* denotes the total proton number density (including 
’bound’ and ’free’ protons). 

In this warm dense matter (WDM) region, correlations and bound state 
formation are relevant. We are considering the influence of the formation of 
light clusters (deuteron d =^H, triton t =^H, helion h =^He, and a =^He) 
on the EOS. Whereas at low densities the nuclear statistical equilibrium 
(NSE) is a reasonable approximation, where non-interacting bound states 
in chemical equilibrium are considered so that the composition is given by 
a mass-action law (continuum contributions are neglected), with increasing 
density the bound states are dissolved due to Pauli blocking and do not 
contribute to the EOS near to the saturation density. 

^Talk at the 32th Max-Born-Symposium and HECOLS workshop on "Three 
Days of Phase Transitions in Compact Stars, Heavy-Ion Collisions and Super¬ 
novae", Institute for Theoretical Physics, University of Wroclaw, Wroclaw, Poland, 
February 17-19, 2014. 
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About thirty years ago a quantum statistical (QS) approach to the 
EOS has been given which allows for a systematic description of correlations 
in nuclear matter. There, the dissolution of bound states because of Pauli 
blocking has been denoted as Mott effect, and the phase diagram shown in 
Fig. 1 has been proposed. Meanwhile, various theoretical investigations have 
been performed as well as hrst experimental evidences [1]. 


Symmetric nuclear matter: Phase diagram 



Figure 1: Phase diagram 
of symmetric nuclear matter 
(schematic), showing the so- 
called Mott line as well as the 
liquid-gas like phase transition 
and the transition to the su¬ 
perfluid state. Presented at 
different talks since 1982, see 
also [5]. 


It is a challenge to nuclear theory to describe the modihcation of the light- 
cluster properties caused by the surrounding nuclear matter. The energies 
Ei(P) of the light clusters to be considered as quasiparticles with energies 
depending not only on the c.o.m. momentum P, but also on the densities 
of neutrons protons and on T. The single nucleon states n, p are 
treated along the same way. 

In principle, one has to start from a many-body Hamiltonian. One has 
to calculate the spectral function for the different channels, and sharp peaks 
can be interpreted as the corresponding quasi-particles. Because we have no 
hrst-principle nucleon-nucleon forces such as the Coulomb force, the choice 
of the interaction potential in the nuclear Hamiltonian remains a problem. 
Different parametrizations are obtained analyzing two-nucleon properties. 
We discuss some simple approximations in the following Section. However, 
to describe also three or four-nucleon systems, the pair potential is no longer 
sufficient, and, for instance, density dependent forces are introduced. 

Therefore it is advantageous to use directly the measured properties such 
as binding energies, rms radii, and scattering phase shifts, avoiding the in¬ 
troduction of a potential. This way, the NSE can be formulated, but also the 
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second virial coefficient can be expressed by measured scattering phase shifts. 
The properties of nuclear matter near saturation density are well reproduced 
by effective energy density functionals such as the Skyrme force or the rela¬ 
tivistic mean-held (RMF) approach with empirical parameter values. These 
single-nucleon quasiparticle energies contain already correlations beyond the 
Hartree-Fock (mean held) approximation. 

To reproduce the EOS, quasiparticle energies for the light clusters can be 
introduced which depend on the thermodynamic variables T, Such 

expressions have been introduced in Refs. PE] and are comprehensively 
discussed in |H]. They can be related to microscopic approaches for the single¬ 
nucleon self-energy and Pauli blocking discussed below. As an alternative, 
the excluded volume concept is also used to mock the Pauli blocking term, 
see P and references given therein. This empirical approach is not directly 
connected with the microscopic QS approach. 


2 Bound-state quasiparticle shifts 

Within a quantum statistical (QS) approach, the EOS ?7,*°*(T,/i„,/Xp) is ex¬ 
pressed in terms of the single-nucleon spectral function which is related to 
the self-energy, see [Tnl ITT] . Performing a cluster decomposition of the self¬ 
energy, the relations 

^ A,u,P 

np°*(T,/i„,/Xp) = ^ , (1) 

^ A,p,P 

are obtained, where P denotes the center of mass (c.o.m.) momentum of the 
cluster (or, for A = 1, the momentum of the nucleon). The internal quantum 
state u contains the proton number Z and neutron number N = A — Z of 
the cluster. 


/AjZ(^) P) hnj hp) r/ AT /y \ /T'l t 1 \A 

exp[(a; - A^/i„ - Zij.p)/T\ - (-1)^ 

is the Bose or Fermi distribution function for even or odd A, respectively, 
that is depending on {T,/r„,/ip}. 
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( 3 ) 


For the A-nucleon cluster, the in-medium Schrodinger equation 

[Er^{pi) H-h Er^{pA) - EA,uiP)]'4’Aup{l ■■■A) 

+ ~ n ^kk''ll^Aup(,l' ... A') = 0 

V...A' i<j kj^i,j 

is derived from the Green function approach. This equation contains the 
effects of the medium in the single-nucleon quasiparticle shift which is given, 
for instance, by the RMF expression |T2j {ns = + n}°^, Yp = n^p^/ ub) 

Er{p; T, ub, Yp) = ^[nirC^ - S'(T, ub, Yp)f + h?c^p‘^ + K(T, hb, Yp) - m^c^. 

(4) 

Approximation formulae for the RMF potentials S, V are found, e.g., in na. 

As well, the effects of the medium are obtained from the Pauli blocking 
terms given by the occupation numbers n(l; T,/i„, Pp) in the phase space 
of single-nucleon states |1) = |pi,(Ti,ri). The occupation numbers can be 
approximated by a Fermi distribution with effective parameter values for 
temperature and chemical potentials, see na. Thus, two effects have to be 
considered, the quasiparticle energy shift and the Pauli blocking. 

We obtain the cluster quasiparticle shifts 

EaAP) - E’iAP) = abIKp) + (5) 

with the free contribution E\^{P) = E\^ + f^P‘^/{2Am). Expressions 
for the in-medium self-energy shift AE^j^{P]T,nB,Yp) and Pauli blocking 
Teff, Ub, are given in [13]. The Coulomb shifts for the light 
elements with Z < 2 considered here are small compared with the other 
contributions and are omitted. 

The Pauli blocking contains the nucleon-nucleon interaction potential 
V{ij, i'j') which has been taken in different approximations. The main fea¬ 
ture is the overlap of the cluster wave function in momentum space with 
the Fermi sphere so that an approximate form of the wave function can be 
used, reproducing characteristic parameters such as the rms radii. Simple 
expressions for the interaction potential have been considered in mm and 
will be discussed in the following Section. 

3 Interaction potentials and Mott line 

As a first step, we consider only A < 2. The generalized Beth-Uhlenbeck 
approach has been investigated by Schmidt et al. n. Results for the sepa¬ 
rable Paris (PEST) interaction potential, rank = 4, are given. The result for 


4 




Figure 2: Two-particle correlations: Frac¬ 
tion of correlations as function of the baryon 
density n for T = 10. Results without in¬ 
medium corrections (dashed) are compared 
with results obtained from the generalized 
Beth-Uhlenbeck formula. Bound state con¬ 
tribution and the continuum contribution 
are also given. From na. 


the composition of symmetric matter at T = 10 MeV is shown in Fig. 2. The 
sharp kink in the bound state contribution to the density at baryon density 
^^Schmidt = 0.013 fm“^ indicates the Mott line where the bound states at 
P = 0 merge with the continuum. The total correlated density Ucorr (con¬ 
taining also the contribution of the continuum) is increasing up to densities 
of about 0.02 fm“^ and then decreasing, vanishing at us 0.1 fm“^. 

Calculations for a more extended parameter range in the T — ub plane 
have been performed by Stein et al. [15]. For simplicity a rank = 1 Yamaguchi 
potential was used which is separable and attractive only. It takes into 
account F-wave scattering (c =^S'o,^5i) and depends on the relative momenta 
of the incoming and outgoing two particles and the coupling strength in the 
respective channel, 


= -A. 




( 6 ) 


where (3 = 1.4488 fm~^ is the inverse potential range, Ai^^, = 2994 MeV fm“^ 
and A35^ = 4264 MeV fm“^ is the coupling strength in the spin-singlet and in 
the triplet channel, respectively. The parameters are fitted to the empirical 
nucleon-nucleon scattering phase shifts and the vacuum bound state energy 
of the deuteron {Ed = —2.225 MeV) which occurs in the spin-triplet channel. 
The coupling of the ^Si to the ^Di channel is neglected. 

Results are shown in Fig. 3. For T = 10 MeV, the Mott line has the 
value n^steini'^ ~ 10 MeV) = 0.02 fm“^ and the maximum contribution of 
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correlated density occurs at 0.04 fm“^, but a contribution of about 20 % 
remains at saturation density. 

We performed calculations within the generalized Beth-Uhlenbeck ap¬ 
proach na for a simple separable Gaussian potential, 




(Pl-P2)^ (Pi-P2)^ 

AcC C 


(7) 


with Arf = 1287.37 MeV for the deuteron (isospin 0) channel, 7 = 1.474 
fm“^, see [101 m, adapted to binding energy and point rms radius of the 
deuteron. The Pauli blocking term for this potential has been evaluated and 
parametrized in HU. For instance, for T = 10 MeV the hrst order Pauli 
blocking shift is jlO] 352.5 MeV fm^ x ub so that the Mott line = 

10 MeV) = 0.0063 fm“^ follows. The exact solution gives 7" = 10 MeV) 
0.009585 fm-T 


4 Different nucleon-nucleon potentials 


To investigate the role of the hard-core repulsion in the nucleon-nucleon 
interaction, we considered a force of the Mongan type |16j 


(Pl-P2)^ (Pl-P2)^ 

V“(12,T2') =-A,,,e ^71— e ^71“ + 


. (Pl-P2)^ 

FT 


(pVp2)^ 
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( 8 ) 

with Ad,a = 1645.89 MeV and ja = 1.749 fm ^ for the attractive part and 
Ad,r = 445.843 MeV and 7 ^ = 2.49 fni“^ for the repulsive part. Only a small 
decrease = 10 MeV) = 0.0094967 fm“^ was obtained. 

The use of a Yamaguchi nucleon-nucleon force with Lorentzian form fac¬ 
tors (j6) gives nearly the same result for the Mott line, l( 7" = 10 MeV) = 
0.0095542 fm-T 


We conclude that Pauli blocking and the Mott line 
are not very sensitive to the details of the nucleon-nucleon interaction poten¬ 
tial. Global features such as the binding energy, the rms radius of the bound 
states as well as the scattering length and the low-/c expansion of the phase 
shifts are characteristics for the extension of the wave function in momen¬ 
tum space. The overlap with the single-particle distribution of the nuclear 
medium is relevant for the Pauli blocking effect. 

The main reason for the different results obtained above for the Mott line 
(’^^Schmidt) ’^^stein) * 1 ^® 1 ° treatment of correlations evaluating the self¬ 
energy and Pauli blocking terms. If in Eq. (|^ the single-nucleon occupation 
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n{i) is replaced by the single-nucleon distribution fi^r[ET{pi,T, ^n, fJ^p)] with 
the chemical potentials and temperature consistent with Eq. Q, only the 
free nucleons are taken into account for the phase occupation. The total 
baryon density ub is larger than the free nucleon density because of the 
contribution of the correlated density. 

It is necessary to take the contribution of correlations into account when 
calculating the mean-field and Pauli blocking shifts. The occupation number 
distribution n{i) can be approximated by a Fermi distribution with parameter 
values describing the phase space occupation by free as well as 

correlated nucleons na. The Mott line obtained this way is = 


Figure 3: Phase diagram of symmetric 
matter, with only two-particle correla¬ 
tions. Fraction of correlations (bound 
state and continuum) to the density are 
given, as well as the Mott line and the 
transition temperatur for superfluid¬ 
ity. From Hg. 


lOMeV) = 0.0087987 fm"’. 


Mon line 



5 Effective medium correction terms 


We give two examples for effective descriptions of the medium modification 
of the contributions of correlations. We can introduce effective binding en¬ 
ergies depending only on the temperature T and the total 

neutron/proton densities nj)’*, , not on the c.o.m. momentum P, and in¬ 

cluding excited states and the continuum of the channel i, using the definition 


Py /^n, Pp)] 

u,P 


= 9 '. 


J2fA.zl-Bf{T,n. 


tot I^tot 




In the non-degenerate case, the sum over P is easily performed, gi is the 
degeneration factor of the ground state. 

Similar to the RMF expressions for the single-nucleon quasiparticle shifts, 
after an appropriate parametrization of np°*) a simple evaluation 

of the EOS ^ is possible. The Mott point = 0 is a signature 
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of medium effects and has been observed recently from measured yields in 
HIC jl]. Theoretical approximations P, 0 E] 77 ,^°** = IQMeV) = 

0.005225 fm“^, = 5MeV) = 0.00292 fm“^ are in reasonable 

agreement with the experimental data |1] = 4.5 MeV) ~ 0.001 

fm“^. 

The other semi-empirical approach is the introduction of an excluded vol¬ 
ume, see p. This parametrization of the medium effects is not very specific 
for the various constituents of the nuclear matter and demands some improve¬ 
ments when it shell be related to the microscopic QS approach. However, it 
is rather efficient for exploratory calculations. 

6 Conclusions 

The quasiparticle energies Ei{P]T,nl^^,n^p^) of the element i depends on 
the c.o.m. momentum P, but also on the thermodynamic parameter values 
for Pauli blocking leads to the dissolution of bound states. In 

particular, = 0 gives the so-called Mott line TMott(^B) The 

correlated part of the density, however, contains besides the contribution of 
bound also the contribution of scattering states. Altogether, the correlated 
density nj(T, np°*) is a smooth function of T, n^°^. 

Calculations for the composition of nuclear matter have been performed 
for the contribution of two-nucleon correlations. We found no strong depen¬ 
dence on the chosen nucleon-nucleon potential. Neglecting the details, one 
can introduce an effective shift of the quasiparticle energies (effective binding 
energy which reproduces the contribution of two-nucleon correlations. 
This is a simple concept to introduce the effect of cluster formation into the 
EOS. For this effective shift, also a Mott point can be introduced where the 
effective shift compensates the binding energy of the cluster. 

Experimental values have been obtained recently [1] and have been com¬ 
pared with theories that use effective shifts EISIH]. Good coincidence has 
been found. However, one should not identify the Mott point obtained there 
with the Mott line describing the dissolution of bound states at P = 0. 

The dependence on temperature cannot be described by a simple excluded 
volume concept. According to IISI. the deuteron Mott line n^°**(T) for sym¬ 
metric matter derived from the first order Pauli blocking shift has the fol¬ 
lowing values: MeV) = 0.0063 fm“^, n^)\“(20MeV) = 0.0128 fm“^, 

(50 MeV) = 0.0385 fm-^, n^(°*|(100 MeV) = 0.0968 fm-^, n^g(140 MeV) 
0.154 fm“^. At very high temperatures deuterons can survive up to satura- 
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tion densities, as also observed recently at LHC experiments. 

Similar approaches are possible also for other clusters with A = 3,4 |13) . 
Because the clusters are more bound, the continuum contributions are less 
relevant. The use of the effective shift instead of the momentum-dependent 
quasiparticle shift (Mott line) is better founded than in the deuteron case. 
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